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Abstract 

The paper discusses local porosity theory and its relation with other 
geometric characterization methods for porous media such as correlation 
functions and contact distributions. Special emphasis is placed on the 
charcterization of geometric observables through Hadwigers theorem in 
stochastic geometry. The four basic Minkowski functionals are introduced 
into local porosity theory, and for the first time a relationship is estab- 
lished between the Euler characteristic and the local percolation prob- 
abilities. Local porosity distributions and local percolation probabilities 
provide a scale dependent characterization of the microstructure of porous 
media that can be used in an effective medium approach to predict trans- 
port. 



1 Introduction 



A crucial prerequisite for the prediction of transport parameters of porous me- 



diais a suitable characterization of the microstructure ( Dulhen, 1992 ; Adler 
1992 ; Sahimi, 1995| ; Hilfer, 199(j ). Despite a long history of scientific study 



the microstructure of porous media continues to be investigated in many ar- 
eas of fundamental and applied research ranging from geophysics ( Hearst ancj 
Nelson, 1985), hydrology (Marsily, 1986; Bear and Verruijt, 1987), petrophysics 
(Lake, 1989) and civil engineering ( Ehlers, 1995|; Diebels and Ehlers, 1996) to 
the materials science of composites (Crivelli-Visconti(ed.), 1998). 

My primary objective in this article is to review briefly the application of 
local porosity theory, introduced in (Hilfer, 1991; Hilfer, 1992; Hilfer, 1996), as 
a method that provides a scale dependent geometric characterization of porous 



or heterogeneous media. A functional theorem of Hadwiger (Hadwiger, 1955, 
p. 39) emphasizes the importance of four set-theoretic functionals for the geomet- 
ric characterization of porous media. In contrast herewith local porosity theory 
has emphasized geometric observables, that are not covered by Hadwigers theo- 
rem ([Hilfer, 1993|; [Hilfer et al., 1994|; [Haslund et al., 1994[). Other theories have 
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stressed the importance of correlation functions (Torquato and Stell, 1982; Stell 



1985) or contact distributions (Levitz and Tchoubar, 1992; Muche and Stoyan, 
1992; 3toyan et al., 1995) for characterization purposes. Recently advances 



in computer and imaging technology have made threedimensional microtomo- 
graphic images more readily available. Exact microscopic solutions are thereby 



becoming possible and have recently been calculated ( 


Widjajakusuma et al. 


1999a; Widjajakusuma et al., 1999b 




Biswal et al., 1999 


). Moreover, the avail- 



ability of threedimensional microstructures allows to test approximate theories 
and geometric models and to distinguish them quantitatively. 

Distinguishing porous microstructures in a quantitative fashion is important 
for reliable predictions and it requires apt geometric observables. Examples 
of im portant geometric observa bles are porosity and specific internal surface 
area ( [Bear, 1972 ; DuUien, 1992 ). It is clear however, that porosity and specific 
internal surface area alone are not sufficient to distinguish the infinite variety 
of porous microstructures. 

Geometrical models for porous media may be roughly subdivided into the 



classical capillary tube and slit models (DuUien, 1992), grain models (Stoyan 
et al., 1995), network models ( Fatt, 1956|; Bryant et al., 1993), percolation mod- 



els (Chatzis and DuUien, 1977; Sahimi, 1993), fractal models (Katz and Thomp- 



son, 1986; Roy and Tarafdar, 1997), stochastic reconstruction models (Quiblier 



1984; Adler, 1992) and diagenctic models (Roberts and Schwartz, 1985; Bakke 



and 0ren, 1997). Little attention is usually paid to match the geometric charac- 



teristics of a model geometry to those of the experimental sample, as witnessed 
by the undiminished popularity of capillary tube models. Usually the match- 
ing of geometric observables is limited to the porosity alone. Recently the idea 
of stochastic reconstruction models has found renewed interest (Adler, 1992; 



Roberts, 1997; Yeong and Torquato, 1998). In stochastic reconstruction models 



one tries to match not only the porosity but also other geometric quantities such 
as specific internal surface, correlation functions, or linear and spherical contact 
distributions. As the number of matched quantities increases one expects that 
also the model approximates better the given sample. Matched models for sed- 
imentary rocks have recently been subjected to a quantitative comparison with 
the experimentally obtained microstructures (Biswal et al., 199£). 



2 Geometrical Problems in Porous Media 

A two-component porous sample § = PUM is defined as the union of two closed 
subsets P C and M C where P denotes the pore space (or component 
1 in a heterogeneous medium) and M denotes the matrix space (or component 
2). For simplicity only two-component media will be considered throughout this 
paper, but most concepts can be generalized to media with an arbitrary finite 
number of components. A particular pore space configuration may be described 
using the characteristic (or indicator) function Xp('") of a set P. It is defined for 
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arbitrary sets P as 



(1 forreP 

[0 forr^P. ^' 

The geometrical problems in porous media arise because in practice the pore 
space configuration Xp('") is usually not known in detail. On the other hand the 
solution of a physical boundary value problem would require detailed knowledge 
of the internal boundary, and hence of X-pi^)- 

While it is becoming feasible to digitize samples of several mm'^ with a 
resolution of a few /im this is not possible for larger samples. For this reason 
the true pore space P is often replaced by a geometric model P. One then 
solves the problem for the model geometry and hopes that its solution u obeys 
M « u in some sense. Such an approach requires quantitative methods for the 
comparison of P and the model P. This in turn raises the problem of finding 
generally applicable quantitative geometric characterization methods that allow 
to evaluate the accuracy of geometric models for porous microstructues. The 
problem of quantitative geometric characterization arises also when one asks 
which geometrical characteristics of the microsctructure P have the greatest 
influence on the properties of the solution u of a given boundary value problem. 

Some authors introduce more than one geometrical model for one and the 
same microstructure when calculating diff'erent physical properties (e.g. diffu- 
sion and conduction). It should be clear that such models make it difficult to 
extract reliable physical or geometrical information. 



3 Geometric Characterizations 
3.1 General Considerations 

A general geometric characterization of stochastic media should provide macro- 
scopic geometric observables that allow to distinguish media with different mi- 
crostructures quantitatively. In general, a stochastic medium is defined as a 
probability distribution on a space of geometries or configurations. Probability 
distributions and expectation values of geometric observables are candidates for 
a general geometric characterization. 

A general geometric characterization should fulfill four criteria to be useful in 



applications. These four criteria were advanced in ( |Hilfcr, 19961) . First, it must 
be well defined. This obvious requirement is sometimes violated. The so called 
"pore size distributions" measured in mercury porosimetry are not geometri- 
cal observables in the sense that they cannot be determined from knowledge of 
the geometry alone. Instead they are capillary pressure curves whose calcula- 
tion involves physical quantities such as surface tension, viscosity or flooding 



history ( [Hilfer, 199q ). Second, the geometric characterization should be di- 
rectly accessible in experiments. The experiments should be independent of 
the quantities to be predicted. Thirdly, the numerical implementation should 
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not require excessive amounts of data. This means that the amount of data 
should be manageable by contemporary data processing technology. Finally, a 
useful geometric characterization should be helpful in the exact or approximate 
theoretical calculations. 



3.2 Geometric Observables 

Well defined geometric observables are the basis for the geometric characteri- 
zation of porous media. A perennial problem in all applications is to identify 
those macroscopic geometric observables that are relevant for distinguishing be- 
tween classes of microstructures. One is interested in those properties of the 
microstructure that influence the macroscopic physical behaviour. In general 
this depends on the details of the physical problem, but some general properties 
of the microstructure such as volume fraction or porosity are known to be rele- 



vant in many situations. Hadwigers theorem (Hadwiger, 1955) is an example of 
a mathematical result that helps to identify an important class of such general 
geometric properties of porous media. It will be seen later, however, that there 
exist important geometric properties that are not members of this class. 

A geometric observable / is a mapping (functional) that assigns to each 
admissible pore space P a real number /(P) — /(P n §) that can be calculated 
from P without solving a physical boundary value problem. A functional whose 
evaluation requires the solution of a physical boundary value problem will be 
called a physical observable. 

Before discussing examples for geometric observables it is necessary to spec- 
ify the admissible geometries P. The set of admissible P is defined as the set 



of all finite unions of compact convex sets ([Hadwiger, 1955 ; Stoyan et al., 1995 



Schneider and Weil, 1992| ; Schneider, 1993 ). Because 3? is closed under unions 



and intersections it is called the convex ring. The choice of is convenient 
for applications because digitized porous media can be considered as elements 
from 51 and because continuous observables defined for convex compact sets can 
be continued to all of Jl. The set of all compact and convex subsets of M.'^ is 
denoted as 3C. For subsequent discussions the Minkowski addition of two sets 
A, B C M'^ is defined as 

A + B = {cc + y : a; e A,y el}. (2) 

Multiplication of A with a scalar is defined by aA = {ax : x E A} for a G M. 

Examples of geometric observables are the volume of P or the surface area 
of the internal dV = dM = FnM. [| Let 

VdiK) = f xJr)d''r (3) 



^The boundary dG of a set G is defined as the difi'erence between the closure and the 
interior of G where the closure is the intersection of all closed sets containing G and the 
interior is the union of all open sets contained in G. 
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denote the d-dimensional Lebesgue volume of the compact convex set K. The 
volume is hence a functional VJ; : IK ^ R on 3C. An example of a compact 
convex set is the unit ball B'' = {a;eR^:la;l<l}=B''(0,l) centered at the 
origin whose volume is 

Other functionals on % can be constructed from the volume by virtue of the 
following fact. For every compact convex K G 3C and every e > there are 
numbers Vj{K),j — 0, . . . ,d depending only on K such that 

d 

Vd{K + el^) = Vj {^)e''-'^id-j (5) 



is a polynomial in s. This result is known as Steiners formula (Hadwiger, 1955 



3toyan et al., 1995). The numbers Vj(K),j = 0. . . ,d define functionals on X 
similar to the volume Vd(]K). The quantities 

wm - ^^^^ (6) 

are called quermassintegrals ( ^chneidcr, 1993| ). From one sees that 



Imi ^{Vd{K + els'*) - Vd{K)) = KiVd-iiK), (7) 

and from that ki = 2. Hence Vd-i(K) may be viewed as half the surface 
area. The functional Vi(K) is related to the mean width w(K) defined as the 
mean value of the distance between a pair of parallel support planes of K. The 
relation is 

V,iK) = ^wiK) (8) 

which reduces to Vi(]K) = w{K)/2 for d = 3. Finally the functional Vb(K) is 
evaluated from (||) by dividing with e'^ and taking the limit e ^ oo. It follows 
that Vb(K) = 1 for all K e 3C \ {0}. One extends Vq to all of % by defining 
Voi^) = 0. The geometric observable Vq is called Euler characteristic. 

The geometric observables Vi have several important properties. They are 
Euclidean invariant (i.e. invariant under rigid motions), additive and monotone. 
Let Td = (M'', +) denote the group of translations with vector addition as group 



operation and let SO{d) be the matrix group of rotations in d dimensions (Barut 



and Raczka, 1986 ). The semidirect product Ed — TdQ SO{d) is the Euclidean 
group of rigid motions in . It is defined as the set of pairs (a, A) with a e T^ 
and A £ SO{d) and group operation 

{a,A)o{b,B)^{a + Ab,AB). (9) 
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An observable / : 3C ^ M is called euclidean invariant or invariant under rigid 
motions if 

fia + AK) = f{K) (10) 

holds for all (a, A) £ Ed and all K e 3C. Here AK = {Ax : x eK} denotes the 
rotation of K and a + K = {a} + K its translation. A geometric observable / is 
called additive if 

/(0) = (11) 
/(Ki U Ka) + /(Ki n K2) = /(Ki) + /(K2) (12) 

holds for all Ki,K2 G OC with Ki U IK2 £ 3C. Finally a functional is called 
monotone if for Ki,]K2 £ % with Ki C IK2 follows /(Ki) < /(K2). 

The special importance of the functionals Vi(K) arises from the following 
theorem of Hadwiger ( Hadwiger, 1955| ). A functional / : 3C ^ K is euclidean 



invariant, additive and monotone if and only if it is a linear combination 

d 



i=0 



with nonnegative constants cq, . . . ,Cd. The condition of monotonicity can be 
replaced with continuity at the expense of allowing also negative and the 
theorem remains vahd ( [Hadwiger, 1955| ). If / is continuous on 3C, additive and 



euclidean invariant it can be additively extended to the convex ring 31 ( schncidci 



and Weil, 1992). The additive extension is unique and given by the inclusion- 



exclusion formula 

/fu^i) = E (-i)i"i-vfnK<) (14) 

\i=l / le9{m) Kiel / 

where J'(m) denotes the family of nonempty subsets of {1, . . . , m} and 1 1 1 is the 
number of elements of I £ J'(m). In particular, the functionals Vi have a unique 



additive extension to the convex ring !K( Schneider and Weil, 1992), which is 
again be denoted by Vi. 

For a threedimensional porous sample with P £ IK the extended functionals 
Vi lead to two frequently used geometric observables. The first is the porosity 
of a porous sample § defined as 

^(Pn§)^03(Pn§) = -^^^^, (15) 

and the second its specific internal surface area which may be defined in view 
of (0) as 
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The two remaining observables (/)i(P n S) Vi(P n S)/y3(S) and n §) = 
Vb(PnS)/V3(§) have received less attention in the porous media literature. The 
Euler characteristic Vq on Ji coincides with the identically named topological 
invariant. For d = 2 and G G one has Vb(G) = c(G) — c'(G) where c(G) is 
the number of connectedness components of G, and c'(G) denotes the number 
of holes (i.e. bounded connectedness components of the complement). 



3.3 Definition of Stochastic Porous Media 

For theoretical purposes the pore space P is frequently viewed as a random set 
( ^toyan et al., 1995 ; Hilfer, 1996| ). In practical applications the pore space is 



usually discretized because of measurement limitations and finite resolution. For 
the purpose of discussion the set S C M'^ is a rectangular parallelepiped whose 
sidelengths are Mi,M2 and M3 in units of the lattice constant a (resolution) 
of a simple cubic lattice. The position vectors = ri-^,,,i^ = (aii,... , ai^) 
with integers 1 < ij < Mj are used to label the lattice points, and is a 
shorthand notation for ri-^,,,i^. Let Yi denote a cubic volume element (voxel) 
centered at the lattice site r^. Then the discretized sample may be represented 
as S = [JiLi "^i- The discretized pore space P, defined as 



u 



V,:, (17) 

{»:X^(rO = l} 

is an approximation to the true pore space P. For simplicity it will be assumed 
that the discretization does not introduce errors, i.e. that P = P, and that each 
voxel is either fully pore or fully matrix. This assumption may be relaxed to 
allow voxel attributes such as internal surface or other quermassintegral den- 
sities. The discretization into voxels reflects the limitations arising from the 
experimental resolution of the porous structure. A discretized pore space for 
a bounded sample belongs to the convex ring 3^ if the voxels are convex and 
compact. Hence, for a simple cubic discretization the pore space belongs to the 
convex ring. A configuration (or microstructure) .Z of a 2-component medium 
may be represented in the simplest case by a sequence 

Z - (Zi, . . . , Z^) = (Xp(ri), . . . , XpivN)) (18) 

where runs through the lattice points and N = M1M2M3. This represen- 
tation corresponds to the simplest discretization in which there are only two 
states for each voxel indicating whether it belongs to pore space or not. In gen- 
eral a voxel could be characterized by more states reflecting the microsctructure 
within the region V^. In the simplest case there is a one-to-one correspondence 
between P and Z given by (|l^). Geometric observables /(P) then correspond 
to functions f{Z) = f{zi, . . . , zn)- 

As a convenient theoretical idealization it is frequently assumed that porous 
media are random realizations drawn from an underlying statistical ensemble. A 
discretized stochastic porous medium is deflned through the discrete probability 
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density 

p{zi, ... ,zn)^ Prob{(Zi = zi) A . . . A (Zw = zn)} (19) 

where Zi € {0, 1} in the simplest case. It should be emphasized that the prob- 
ability density p is mainly of theoretical interest. In practice p is usually not 
known. An infinitely extended medium or microstructure is called stationary or 
statistically homogeneous if p is invariant under spatial translations. It is called 
isotropic Hp is invariant under rotations. 

3.4 Moment Functions and Correlation Functions 

A stochastic medium was defined through its probability distribution p. In prac- 
tice p will be even less accessible than the microstructure ¥ = Z itself. Partial 
information about p can be obtained by measuring or calculating expectation 
values of a geometric observable /. These are defined as 

1 1 

{f{zi, ... , zn)) = X! ■ • ■ X! ■^'^^1' ■ ■ ■ ' ^n)p{zi, ... ,zn) (20) 

2:1=0 zn—O 

where the summations indicate a summation over all configurations. Consider 
for example the porosity (/)(S) defined in For a stochastic medium 0(§) 

becomes a random variable. Its expectation is 

1/3(8) Vsin) J '/ 

s 

N N 

^ ' Z— 1 2 — 1 

1 ^ 1 ^ 

= ]^ E = 1} = ]^ E e P} 

i=\ 1=1 

If the medium is statistically homogeneous then 

(0) = Prob{z, = 1} = Prob{r, e P} = (xp(n)) (22) 

independent oil. It happens frequently that one is given only a single sample, 
not an ensemble of samples. It is then necessary to invoke an ergodic hypothesis 
that allows to equate spatial averages with ensemble averages. 

The porosity is the first member in a hierarchy of moment functions. The 
n-th order moment function is defined generally as 

5„(ri,... ,r„) = (xp(r-i)...Xp(r-«)) (23) 



8 



(21) 



for n < TV. For stationary media 5'„(ri, . . . r„) = g{ri — Vn, ■ • ■ , r„_i — r„) 
where the function g depends only on 71 — 1 variables. Another frequently 
used expectation value is the correlation function which is related to 82- For a 
homogeneous medium it is defined as 

G(.o, r) . Gir - ^ i^^^^^^^M ~ i^f _ S.ir - .0) - (^r (ro))^ 



-(</>)) ^i(ro)(l-^i(ro)) 

(24) 

where rg is an arbitrary reference point, and (0) — S'i(ro). If the medium is 
isotropic then G(r) = G(|r|) = G{r). Note that G is normalized such that 
G(0) = 1 and G(oo) = 0. 

The hierarchy of moment functions S„, similar to p, is mainly of theoretical 
interest. For a homogeneous medium S„ is a function of n — 1 variables. To 
specify Sn numerically becomes impractical as n increases. If only 100 points 
are required along each coordinate axis then giving Sn would require IQ^''^""!) 
numbers. For d = 3 this implies that already at n = 3 it becomes economical to 
specify the microstructure P directly rather than incompletely through moment 
or correlation functions. 



3.5 Contact Distributions 



An interesting geometric characteristic introduced and discussed in the field 



of stochastic geometry are contact distributions (Delfiner, 1972; Stoyan et al 
206). Certain special cases of contact distributions have appeared 
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also in the porous media literature (DuUicn, 1992). Let G be a compact test 
set containing the origin 0. Then the contact distribution is defined as the 
conditional probability 



HG{r) = 1 - Prob{0 ^ M + (- 



I} = 1 



Prob{M n rG = 0} 



(25) 



If one defines the random variable R = inf{s : M n sG ^0} then HQ{r) = 
Prob{i? < r\R > 0} ( [Stoyan et al., 1995] ). 

For the unit ball G = B(0, 1) in three dimensions Hm is called spherical 
contact distribution. The quantity 1 ~ H^{r) is then the distribution function of 
the random distance from a randomly chosen point in P to its nearest neighbour 
in M. The probability density 



ar ar 



(26) 



was discussed in ( ^cheidegger, 1974 ) as a well defined alternative to the fre- 
quently used pore size distrubution from mercury porosimetry. 

If a voxel has other attributes besides being pore or matrix one may define also mixed mo- 
ment functions Si-^...i„ (ri, . . . , r„) = (ri) (r„)) where 0i{rj) = Vi{FnYj)/Vi{Yj) 
for i = 1 , . . . d are the quermassintegral densities for the voxel at site r j . 
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For an oriented unit interval G = B^(0, 1; e) where e is the unit vector one 
obtains the hnear contact distribution. The linear contact distribution written 
as L{re) = 0(1 — ^^Bi(o.i:e) ('')) is sometimes called lineal path function (Ycong 



and Torquato, 1998). It is related to the chord length distribution pc;(^) defined 
as the probability that an interval in the intersection of P with a straight line 
containing B^(0, l;e) has length smaller than x (Hilfcr, 199£; [Stoyan et al. 



1995, p. 208). 



3.6 Local Porosity Distributions 

The idea of local porosity distributions is to measure geometric observables 
inside compact convex subsets IK C §, and to collect the results into empirical 



histograms (Hilfer. 1991). Let IK(r, L) denote a cube of side length L centered at 



the lattice vector r. The set IK(r,L) is called a measurement cell. A geometric 
observable /, when measured inside a measurement cell K(r, L), is denoted as 
/(r, L) and called a local observable. An example are local Hadwiger functional 
densities / = J2i=Q ^i'^i with coefficients q as in Hadwigers theorem (|l^). Here 
the local quermassintegrals are defined using as 



V',(PnlK(r,i)) 



VK,(PnK(r,i)) 
Vd(K(r,L)) 



for i = 0, . . . , d. In the following mainly the special case d 
For = 3 the local porosity is defined by setting i = 0, 

0(r,L) = Vo(Fn]K(r,L)). 



(27) 



3 will be of interest. 



(28) 



Local densities of surface area, mean curvature and Euler characteristic may be 
defined analogously. The local porosity distribution, defined as 



M0;r,i) = (5(0-0(r,L))) 



(29) 



gives the probability density to find a local porosity 0(r, L) in the measurement 
cell IK(r,L). Here 5{x) denotes the Dirac ^-distribution. The support of ^ 
is the unit interval. For noncubic measurement cells IK one defines analogously 
/x(0; K) = {5{(j) - (j){K))) where (/>(K) = 0(PnlK) is the local observ able in cell K 
The concept of local porosity distributio ns P| was intro duced in ( Hilfcr, 1991 ) 
and has been generalized in two directions (Hilfer, 1996). Firstly by admitting 
more than one measurement cell, and secondly by admitting more than one 
geometric observable. The general n-cell distribution function is defined as 
([Hilfer, 199^) 



Mn;/i,... ,/„ (/1I7 • • • , fln\ ■ ■ ■ ] fnl, ■ ■ ■ , fnrn: Ki, . . . , IK„) — 

{S{f 11 - .fl(IKi)) . . . (5(/i„ - /i(IK„)) . . . 6{Ul - /l(IKi)) . . . S{f,nn - /m(K„))) 

(30) 



or more generally "local geometry distributions" ( Hilfer, 1995 ; Hilfcr, 199t ) 
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for n general measureraent cells Ki, . . . ,K„ and m observables /i, . . . , /,„. The 
n-cell distribution is the probability density to find the values /n of the local 
observable /i in cell Ki and /12 in cell K2 and so on until fmn of local observable 
fm in K„. Definition ( |30| ) is a broad generalization of (p9[). This generalization 
is not purely academic, but was motivated by problems of fluid flow in porous 
media where not only ipo but also ipi becomes important ( Hilfcr, 1992| ). Local 



quermassintegrals, defined in (p7|), and their linear combinations (Hadwiger 
functionals) furnish important examples for local observa bles in (^0|), and they 



have recently been measured on real sandstone samples ( Manwart and Hilfcrj 



19991) 



The general n-cell distribution in (30) is very general indeed. It even contains 



p from (19) as the special case m — 1, fi — (p and n — N with 



K(ri,a). More precisely one has 

^J.N■,pi<t'l, ■ • ■ ,0jv; Vi, . . . ,Yn) = pi<Pi, ■ • ■ An) (31) 

because in that case 0i = = 1 if Vi G P and (pi — Zi — ioi Y ^ P. In 
this way it is seen that the very definition of a stochastic geometry is related 
to local porosity distributions (or more generally local geometry distributions). 
As a consequence the general n-cell distribution /^n;/i,... ,/„ is again mainly of 
theoretical interest, and usually unavailable for practical computations. 

Expectation values with respect to p have generalizations to averages with 
respect to /z. Averaging with respect to /i will be denoted by an overlinc. In 
t he special c ase m = 1, fi =0 and Ki = Vi = K{ri, a) with n < N one finds 
( iHilfer, 199^ ) 



{ri,a) ■ ■ ■ (p[rn,a) 
1 1 




1 1 



;Vi, . . . ,V„)d0i • ■ -dcpn 
4>i ■ ■ ■ 4'n^J■N■,<pi4'l, ■ ■ ■ ,0Ar; Vi, . . . ,VAr)d0i • ■ -dtpN 




1 1 



01 • ■ • </>n (^(01 - 0(t-i, a)) • • • 5{(I)N - 0(rAr, a))) d0i • ■ • d0jv 

= (0('^i,a) • • •0(r„,a)) 

= (Xp(^l)---Xp(^n)) 

= 5„(ri,... ,r„) (32) 

thereby identifying the moment functions of order n as averages with respect to 
an n-cell distribution. 

For practical applications the 1-cell local porosity distributions /i(r, L) and 
their analogues for other quermassintegrals are of greatest interest. For a ho- 
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mogeneous medium the local porosity distribution obeys 



H{(f>; r, L) = ^{(f>; 0, L) = ^{(f>; L) 



(33) 



for all lattice vectors r, i.e. it is independent of the placement of the measure- 
ment cell. A disordered medium with substitutional disorder ( ^iman, 1982| ) may 
be viewed as a stochastic geometry obtained by placing random elements at the 
cells or sites of a fixed regular substitution lattice. For a substitutionally disor- 
dered medium the local porosity distribution /i(r', L) is a periodic function of r 
whose period is the lattice constant of the substitution lattice. For stereological 
issues in the measurement of n from thin sections see (Virgin et al., 1996). 

Averages with respect to /x are denoted by an overline. For a homogeneous 
medium the average local porosity is found as 



1 

(r,L) = j M0;r,L)d0=(<^)=^ 



(34) 



independent of r and L. The variance of local porosities for a homogeneous 
medium defined in the first equality 



L3 



(0) (!-('/'» 1 



2 

l3 



E ^(n-r,) 



(35) 



is related to the correlation function as given in the second equality ( HilferJ 



1996 ). The skewness of the local porosity distribution is defined as the average 

(36) 



a{Lf 



The limits L — > and i — > oo of small and large measurement cells are of 
special interest. In the first case one reaches the limiting resolution at L — a 
and finds for a homogeneous medium ( ]Hilfer, 1991 ; Hilfer, 1996) 



M(0;a)-0<5(0-l)-(l-(/.)(5(0). 



(37) 



The limit L ^ oo is more intricate because it requires also the limit § ^ M"^. 
For a homogeneous medium ( ^5|) shows (t{L) ^ for L ^ and this suggests 

^i{(|3,L ^ oo) = S {(!}-'$). (38) 

For macroscopically heterogeneous media, however, the limiting distribution 



may deviate from this result (Hilfer, 1996). If (38) holds then in both limits the 
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geometrical information contained in /i reduces to the single number (f) = (0). 
If (|3^) and ( p8| ) hold there exists a special length scale L* defined as 

L* = min{L : ^(0; L) = fi{l; L) = 0} (39) 

at which the (5-coniponents at </) = and (p — 1 vanish. The length L* is a 
measure for the size of pores. 

The ensemble picture underlying the definition of a stochastic medium is 
an idealization. In practice one is given only a single realization and has to 
resort to an ergodic hypothesis for obtaining an estimate of the local porosity 
distributions. The local porosity distribution may then be estimated by 

/I(0;i) = lV<5(0-0(r,L)) (40) 

TO ^ — ' 

r 

where m is the number of placements of the measurement cell K(r, L). Ideally 
the measurement cells should be far apart or at least nonoverlapping, but in 
practice this restriction cannot be observed because the samples are not large 
enough. The use of ]! instead of can lead to deviations due to violations of 
t he ergodic hypothesis or simply due t o ovcrsampling the central regions of § 



( Piswal ct al., 1998| ; [Biswal et al., 19991) 



3.7 Local Percolation Probabilities 

Transport and propagation in porous media are controlled by the connectivity 
of the pore space. Local percolation probabilities characterize the connectivity 



(Hilfcr, 1991). Their calculation requires a threedimensional pore space repre- 



sentation, and early results were restricted to samples reconstructed laboriously 



from sequential thin sectioning (Hilfer et al., 1997) In this section a relation 



ship between the Euler characteristic and the local percolation probabililties is 
established for the first time. 

Consider the functional A : 3C x 3C x 31 ^ Z2 = {0, 1} defined by 

f 1 : if Ko Koc^ in P 
I : otherwise 

where IKq C M"^, Koo C M"^ are two compact convex sets with Kon(Pn§) ^ and 
n (P n S) ^ 0, and "Kq ^ Kqo in P" means that there is a path connecting 
Kg and Kqo that lies completely in P. In the examples below the sets Kg and 
Koo correspond to opposite faces of the sample, but in general other choices are 
allowed. Analogous to A, which is defined for the whole sample, one defines for 
a measurement cell 

A„(r, L) = A(Ko„, K^„; P n K(r, L)) = ^ ^ ''^'^ 7" ^ (42) 

I : otherwise 

where a — x,y^ z and Kqj,, '^oox denote those two faces of K(r, L) that are nor- 
mal to the X direction. Similarly IKq^ , Kqoj,, IKq^Kooz denote the faces of K(r, L) 
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normal to the y- and z-directions. Two additional percolation observables A3 
and Ac are introduced by 



A3(r, L) = A,(r, L)Ky{r, L)A,{r, L) (43) 
Ac(r, L) = sgn(A,(r, L) + A^(r, L) + A,(r, L)). (44) 

A3 indicates that the cell is percolating in all three directions while Ac indicates 
percolation in x- or y- or z-direction. The local percolation probabilities are 
defined as 

A„(</.;L) = ^4^%^:lte^ (45) 

where 

fl: if,^ = 0(r,L) 
^^•^('■'^) = \0: otherwise. ^''^ 

The local percolation probability Aq, (0; L) gives the fraction of measurement 
cells of sidelength L with local porosity (j) that are percolating in the "a"- 
direction. The total fraction of cells percolating along the "a" -direction is then 
obtained by integration 



p^{L) = / /x(</.;L)A„(0;L)d0. (47) 



This geometric observable is a quantitative measure for the number of elements 

that have to be percolating if the pore space geometry is approximated by a 
substitutionally disordered lattice or network model. Note that neither A nor 
A„ are additive functionals, and hence local percolation probabilities are not 
covered by Hadwigers theorem. 

It is interesting that there is a relation between the local percolation prob- 
abilities and the local Euler characteristic Vb(P fl lK(r, I)). The relation arises 
from the observation that the voxels are closed, convex sets, and hence for 
any two voxels , Yj the Euler characteristic of their intersection 

Wnv,) = ^^^^^"^^-^! (48) 

^ [0 : if ViOVj =0 ^ ^ 

indicates whether two voxels are nearest neighbours. A measurement cell IK(r, L) 
contains voxels. It is then possible to construct a {L^ -|- 2) x {L^ + 2)^-matrix 
B with matrix elements 

(i3),(,,,) =t/o(V,nV,) (49) 

(B),(,,,) = -yo(v.nv,) (50) 
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where i,j G {0, 1, ... , L^, 00} and the sets Vq = Kg and Voo = ^00 are two 
opposite faces of the measurement ceh. The rows in the matrix B correspond to 
voxels while the columns correspond to voxel pairs. Define the matrix A = BB^ 
where B^ is the transpose of B. The diagonal elements (^4)^^ give the number of 
voxels to which the voxel Vi is connected. A matrix element {A)ij differs from 
zero if and only if Vi and Vj are connected. Hence the matrix A reflects the 
local connectedness of the pore space around a single voxel. Sufficiently high 
powers of A provide information about the global connectedness of P. One finds 

A(Ko,IK^;PnK(r,L)) = sgn(|(A")ooo|) (51) 

where (A'")ooo is the matrix element in the upper right hand corner and m is 
arbitrary subject to the condition m > Li^ . The set P n K(r, L) can always be 
decomposed uniquely into pairwise disjoint connectedness components (clusters) 
Bi whose number is given by the rank of B. Hence 

ranks 

K,(Fn]K(r,L))= (52) 
1=1 

provides an indirect connection between the local Euler characteristic and the 
local percolation probabilities mediated by the matrix B. ^) 

The theoretical concepts for the geometric characterization of porous media 
discussed here are also useful in effective medium calculations of transport pa- 



rameters such as conductivity or permeability (Hilfer, 1991; Hilfcr, 1992; Hilfcr 



1996). The resulting parameterfree predictions agree well with the exact result 



( Widjajakusuma et al., 1999a; Widjajakusuma et al., 1999b). 



4 Application Example 



The purpose of this section is to show that the previous theoretical framework 
can be used directly for predicting transport in porous media quantitatively and 
without free macroscopic fit parameters. The theory presented above allows a 
quantitative micro-macro transition in porous media and opens the possibility 
to determine the elusive "representative elementary volume" needed for macro- 
scopic theories in a quantitative and property specific manner. 



In Reference (Biswal et al., 1999) a fully threedimensional experimental sam- 
ple of Fontainebleau sandstone was compared with three geometric models, some 
of which had not only the same porosity and specific internal surface area but 
also the same correlation function G{r). A large number of the geometric quan- 



tities discussed above was calculated in (Biswal et al., 1999). For a discussion 
on the influence contact of distrtributions see ( Manwart et al., 2000] ). The total 



fraction of percolating cells, defined in eq. ( p^ ) above, was measured in (Biswal 
st al., 1999) for Fontainebleau and some of its models. Figure |l| shows the 



*For percolation systems it has been conjectured that the zero of the Euler characteristic 
as a function of the occup ation probability is an approximation to the percolation threshold 



(Mecke and Wagner, 1991) 
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Figure 1: Total fraction of percolating cells for Fontaine bleau sandstone (EX ) 
and three of its models (DM,GF,SA) as discussed in Ref. (Biswal et al., 1999) 



total fraction of percolating cells as a function of length scale (side length of 
measurement cells). 

It turns out that the quantity (L) displayed in Figure |l| correlates very well 
with transport properties such as the hydraulic permeability. Recently transport 
properties such as the permesbilities and formation factors of the Fontainebleau 
sandstone and its geometries were calculated numerically exactly by solving the 
appropriate microscopic equations of motion on the computer ( Manwart ancj 
Hilfer, 1999| ). Some of the results are summarized in Table |l| below. 

One sees from Table Q that while EX and DM are very similar in their 
permeabilities and formation factors the samples EX and GF have significantly 
lower values with GF being somewhat higher than SA. The same relationship 
is observed in Figure |l] for the percolation properties. These results show that 
the purely geometrical local percolation probabilities correlate surprisingly well 
with hydraulic permeability and electrical conductivity that determine physical 
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Table 1: Physical transport properties of Fontainebleau sandstone and three 
geometric models for it (see ( Biswal et al., 1999D ). an is the conductivity in 
the direction i = x,y,z in units of 10~"^crp, where crp is the conductivity of 
thematerial filling the pore space, ku is the permeability in the direction i = 
X, y, z in mD. 





EX 


DM 


SA 


GF 


kzz [mD] 


692 


923 


35 


34 


kyy[mD] 


911 


581 


22 


35 


k^^[mD] 


790 


623 


20 


36 




18.5 


26.2 


1.35 


2.05 


ayy[lO~^ar] 


21.9 


17.0 


0.87 


1.97 




20.5 


17.1 


0.96 


1.98 



transport. 
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